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ABSTRACT 



The effect of injection of up to ten percent by volume air in damping 
the amplitude of water borne pressure pulses was studied. Pressure pulses 
were generated using a water hammer model. The water hammer model con- 
sisted of a length of pipe with water flowing through it and a quick clos- 
ing valve located downstream. Measured amounts of air were injected into 
a water stream and the valve was quickly closed. The pres sure “time history 
was recorded on oscilloscope photographs at two locations along the pipe. 

A mathematical model was developed for this two phase model. Small amounts 
of air greatly reduced the amplitude of the pressure pulses. By injection 
of one percent by volume of air into the model at 15.5 psia, the pressure 
amplitude was reduced to 8.6 percent of the theoretical no air amplitude. 
The effect decreased with increased pressure of the system. 
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Introduct ion* 



In recent years the United States Navy has been doing a great deal 
of research and development in methods of reducing the noise level of 
machinery for submarines and for anti-submarine ships. The fluid systems 
of these ships transmit noises via pressure pulses in the systems* con- 
tents, A similar problem concerning fluid borne pressure pulses is that 
of vibrations set up by water hammer in piping systems. The objective of 
this study was to quantitatively determine the effect of air in damping 
water borne pressure pulses. 
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2. Theory. 

To study the effect of air in damping water borne pressure pulses, 
it was necessary to generate a pressure pulse in the liquid. For this 
study, the water hammer model was used to generate the pressure pulse. 

The water hammer model shown in Fig. 1 consists essentially of a head 
tank, a length of pipe, and a quick closing valve. The water is allowed 
to flow through the pipe with a velocity head. If the valve is quickly 
closed, this velocity head is immediately changed to pressure head which 
in turn sends a pressure pulse down the pipe. This pressure wave travels 
down the pipe at the speed of sound in the fluid contained in the system. 

Considering the system shown in Fig. 1, the maximum pressure rise by 
instantaneously closing the valve is given by the following equation. 

Where: H = Head rise above the Initial static no flow head, ft. of water 

y 

= Velocity of the water in the pipe before valve closure, ft/sec. 

g = Proportionality constant in Newton’s second law of motion, 

° (lbm-ft)/ (lbf-sec 2 ). 

a = Velocity of the pressure pulse, ft/sec. 

The derivation of this equation as derived by Stepanoff / 1 3/ is included 
in Appendix B. This equation does not include frictional effects. If 
the valve is closed instantaneously, the maximum pressure given by equa- 
tion (1) will be the same along the whole length of the pipe. Obviously, 
a valve cannot be closed instantaneously. However, this same maximum 
pressure rise at a given location A along the pipe will be attained if 
the valve is closed in less time than it takes the pressure pulse to travel 
from position A to the tank and to reflect back to position A. At the 
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valve, the time for the pressure pulse to reflect back is 2h/a. Where 
L is the length of the pipe from the valve to the tank. At the center 
of the pipe, the time for the pressure pulse to reflect back is L/a . 

Next to the tank the maximum pressure will not be attained unless the 
valve is closed instantaneously. For the maximum pressure given by 
equation (1) to be attained at a given point located at a distance X 
from the valve, the closure time must be less than 2.(L-X)/d * 

The speed of the pressure wave in a pure water pipe system is given 
by 

^ 1 

( 2 ) 



a = 



Vf* ( '/kw 4 egr) 

Where: = Water density, slugs/ft^ 

2 

K = Water bulk modulus, Ibf/ft . 
w 

D = Inside pipe diameter, in. 
e = Pipe wall thickness, in. 

E = Pipe modulus of elasticity. 

This equation has been derived in Appendix B. This equation can be ex- 
tended to a two phase air-water system. The speed of a pressure pulse 
based upon an isothermal compression of the air by the pressure pulse is 
given by 



a i “ qj f + X( '-*) \p . + e % [0-^1 



(3) 



Where: x = Volume fraction air at water inlet conditions, 

a^. = Velocity of the pressure pulse in the air-water mixture 
with isothermal compression of the air by the pressure 
pulse, ft/sec. 

a = Velocity of sound in air at the temperature of the fluid, 
a ft/sec. 
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a = Velocity of sound in water at the temperature of the fluid, 



w 



ft/sec. 



k = Ratio of specific heats at constant pressure to that at 
constant volume for air, 

2 

P = Absolute pressure of the system, Ibf/ft . 

ABS 



a = Air density, slugs/ft ^ 

°)^ = Compressibility of water, ft^/lbf. 



The speed of the pressure pulse based upon an adiabatic compression of the 
air by the pressure pulse is given by: 




X 

q 



Q z . 






(3a) 



Equations (3) and (3a) have been derived in Appendix B. It is noted that 
equations (3) and (3a) reduce to equation (2) when no air is present. 
Substituting equation (2) into equation (1) and multiplying by the specific 
weight of water ( Y^) , gives the maximum pressure rise upon valve closure 
with no air in the pipe. 







Vo 



X 



LuL. 



x/psv (/K 



W 




(4) 



The maximum pressure rise for an air-water mixture based upon isothermal 
compression of air is given by: 



P - 




(5) 



The term a^ is given by equation (3). The derivation of equation (5) is 
in Appendix B. The same equation based upon adiabatic compression of the 
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air is given by: 



P = 



Vo T c a 



c 



9 * 



(5a) 



The term a is given by equation (3a). The derivation of equation (5a) 

A 

is in Appendix B. The ratio of equation (5) to equation (4) gives: 




The ratio of equation (5a) to equation (4) gives: 




If equation (6a) is used, adiabatic conditions are assumed or if equation 
(6) is used, isothermal conditions are assumed. Normally, sound traveling 
in the atmosphere is adiabatic. Karplus / 4/ conducted experiments with 
air bubbles approximately 0.01 centimeters in diameter in an air-water 
mixture and found that his experimental data was best correlated by assum- 
ing an isothermal compression. The method used to inject air for this 
investigation did not give the small bubble distribution Karplus obtained. 
In Karplus 1 method, there was a high heat transfer surface per unit bubble 
volume. Therefore, the compression was nearly isothermal in the water. 
With larger bubbles there is less surface per unit bubble volume and the 
situation is probably best described by the adiabatic equation (6a). 

Equation (6a) provides the fundamental theory for the damping of 
pressure pulses in water by air injection. By measuring the air injection 
rate, the water flow rate, and the pressure-time history at a given loca- 
tion in the system of Fig. 1, the effect of air in damping a pressure 
pulse in water is obtained. In this report, the experimental pressures 
obtained are compared with the theoretical relationship of equation (6a). 
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3. Equipment. 

To study the effect of air in damping water borne pressure pulses, 
and to verify the applicability of equation (6a) to the water hammer 
model, the system shown in Fig. 1 was constructed. 

The system consisted of 1) a head tank, a pipe and a weigh tank 
system 2) an air injection system 3) a valve and closure system and 
4) a pressure recording system. A general description for each part of the 
system will be presented here. Exact detailed equipment specifications are 
included in Appendix C. 

The head tank was a tank two feet in diameter with four constant head 
overflow sections. These overflow sections were located nominally at two, 
four, six and eight feet above the main tank outlet line. The main pipe 
length consisted of 59 feet 6 and 7/8 inches of two inch extra strong 
black steel pipe. This pipe was supported at ten feet intervals by a timber 
structure. At the main pipe discharge, a 20 inch by 28 inch by 20 inch 
weigh tank, platform scale and electric timer were used to measure the 
weight flow rate of the water. Water for the system was supplied from the 
city water supply. 

The air injection system contained a Fischer and Porter Company roto- 
meter with a glass ball float for air flow measurement. This rotometer 
had a maximum capacity of 0.325 cubic feet per minute at 14.7 psia and 
70°F. Compressed air from a nominal 100 psig compressor was reduced by a 
regulating valve to approximately fifteen psig, A mercury manometer was 
used to measure the air pressure, and the air temperature was measured 
with a copper constantan thermocouple junction upstream of the rotometer. 
The air was transferred from the regulating valve to the point of injection 
through one quarter inch vinyl tubing. The air was regulated by a manual 
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control valve located immediately following the rotometer. Fig. 1 shows 
a diagram of the system. Fig. 2 shows a photograph of the component parts 
of the air injection system. The air was injected into the pipe near the 
tank through a system of three circles of small holes drilled equally 
spaced around the circumference of the pipe. Details are included in 
Appendix C. A concentric plexiglass cylinder with "0" ring seals surrounds 
these holes. The system was designed so that any one, any two, or any three 
of the circles of holes could be used to inject air at any one time. Fig. 

5 shows a cross sectional view of the air inject annular ring. This method 
for air injection is very simple and could easily be adapted in practical 
applications . 

If the same maximum pressure (neglecting friction losses) was to be 
obtained along one half the length of the pipe, it was necessary to close 
the valve in less than one half the pipe period This was obtained 

by the use of a modified two-inch Crane quick-closing valve. A steel cy- 
linder two inches in diameter and one inch high was screwed onto the top 
of the stem. The valve was closed by dropping a seven pound weight from 
nine feet onto the cylinder. A spring loaded catch was provided to keep 
the valve closed when the valve was in the shut position. An isolated 
stop stand was provided for stopping the valve. This step stand minimized 
the effect of the valve stopping impact from being transmitted to the pipe 
and to the valve body. A quick release mechanism was provided to release 
the weight which was guided down through a perforated aluminum tube, one 
and one-quarter inches in diameter. Fig. 3 shows a photograph of the 
valve and closure system. In Appendix G, calculations are provided which 
were used to estimate the closure time of the valve. The closure time of 
the valve was measured with an oscilloscope and a Polaroid camera. The 



7 



time sweep was triggered with a photoelectric cell. The closure was 
marked by a 45 volt direct current signal on the vertical plates of the 
oscilloscope. The distance that the sweep moved on the oscilloscope 
from trigger to closure, was photographed with 3000 speed Polaroid film. 

The closure time for the seven pound weight from nine feet was 0,007 + 0,003 
seconds. This agrees closely with the values calculated in Appendix G. 

A Kistler model 401 quartz pressure transducer was installed at three 
locations along the pipe, A Kistler 407-5 low noise cable was used to 
connect the transducer to the Kistler model 651 piezo-calibrator amplifier, 

A shielded cable was used to connect the output of the piezo-calibrator 
amplifier to a Hewlett Packard model 130 oscilloscope. An oscilloscope 
Polaroid camera was used with Polaroid 3000 speed film to take a time 
exposure of the pressure- time history for each run. When the valve was 
completely closed, a 45 volt direct current circuit was completed. This 
triggered the oscilloscope time sweep. Fig. 4 is a photograph of the com- 
ponent parts of the pressure recording system. 
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4. 



Procedure . 



The objective of the experimental part of this work was to verify 
the applicability of the theoretical model. 

The pressure-time history was recorded at three locations along the 
pipe. Fig. 13 shows these locations. They were numbered starting from 
the valve, one, two, and three respectively. The pressure-time history 
was recorded at each transducer location for five different air-water 
mixtures at each nominal head in the head tank of two* four, six and 
eight feet. Since the effect of small amounts of air was of interest, 
air rates were chosen such that the volume percent air in the system was 
always less than ten percent. 

Because of the effect of as little as 0.1 percent by volume air in 
the water for the no air runs, it was difficult to get good results with 
water taken directly from the water main. For these runs, it was found 
that better results were obtained if the water was left in the tank over- 
night. 

The transducer was calibrated using an Ashcroft dead-weight tester. 

The calibration curves are given in Appendix C, Fig. 15 and 16. 

A calibration curve for the rotometer was available from the manu- 
facturer. Four points on this curve were checked with a positive displace- 
ment wet gas meter. The points checked closely, and the curve was used. 

The detailed procedure for each run is given in Appendix D and a 
sample data sheet is shown in Appendix E. 

Only one pressure transducer was used. This transducer was moved 
to the desired location for each run. 
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5. 



Results . 



The results have been presented in the form of tables, photographs, 
and curves. In Table 1 ’'transducer locations" are numbered as shown in 
Appendix A, Fig. 13. The "head” refers to the head of water in the head 
tank. The values given are nominal, and the exact values of the heads are 
shown in Appendix A, Fig. 13. The term Y~| is the volume percent of air 
at the water inlet conditions of temperature and pressure. The term 
is the weight percent of air in the mixture. P is the maximum pressure 
rise above the static pressure with no flow in the pipe. The ratio P/P^, 
is the ratio of the maximum pressure obtained when the valve is closed to 
the theoretical maximum pressure with no air, equation (4). 

Results are tabulated for transducer locations one and two. Meaningful 
results were unattainable from location three. The reasons for this will be 
discussed later. 

Fig. 6 shows a series of photographs for nominal heads of two and four 
feet of water and Fig. 7 shows the series of photographs for nominal heads 
of six and eight feet of water. Each row of photographs is for a given 
head. In the first photograph no air was injected into the water. The 
amount of air increases moving from left to right. The camera used to take 
the photographs takes the pictures such that the sweep of the oscilloscope 
moves from right to left rather than left to right as the oscilloscope 
sweep actually moves. The initial point in Fig. 6 and 7 has been marked by 
lines. The horizontal arrows point in the direction that the sweep moves 
along the horizontal scale. The vertical arrows point in the direction of 
increasing pressure. In the photographs, it should be noted that the pres- 
sure and time scales are different when air Is Injected from those when no 
air is injected. 
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In Fig. 8 through 11, the volume percent of air was based upon 
the water conditions of temperature and pressure at the inlet to the main 
run of pipe. Fig. 8, 9, 10, and 11 are respectively for nominal heads of 
two, four, six and eight feet of water. Both the theoretical adiabatic 
(Equation (6a)) and theoretical isothermal (Equation (6)) curves have been 
drawn in. The experimental points for transducer locations one and two are 
plotted and marked with symbols. 



11 



6 . 



Discussion of Results, 



Transducer location three could not be used to obtain meaningful 
data because the valve could not close fast enough. The valve would 
have had to close in less than 0.0015 seconds to get the maximum pressure 

a' 

without interference of the reflected pulse returning from the tank. Hence 
only data from transducer locations one and two are presented. This sub- 
ject was previously discussed in the theory section. 

As the tabulated results show, two runs were made for each set of run 
conditions. This was done as a means of checking the accuracy and consist- 
ency of the results. 

The weight percent of air in the two phase water-air mixture is includ 
ed in the tabulated results to show the small weight fraction of air pre- 
sent in the air-water mixtures. 

It was impossible to obtain meaningful results at transducer location 
one with no air in the water. The high frequency pressure pulses set up 
by the valve stem and gate vibrating made it impossible to distinguish 
between the pulse created by stopping the water and that created by the 
vibrating valve parts. With a small amount of air, these high frequency 
pulses were damped out and meaningful results were obtained. The effect 
of air in damping the high frequency noise is seen in Fig. 6 and 7. 

The pressure-t ime photographs of Fig. 6 and 7 show how the pressure- 
time history at transducer location two appeared on the oscilloscope. The 
photographs show that in the region observed; 1) the greater the percent 
of air present in the mixture, the smaller the pressure amplitude and the 
slower the pressure pulse moves down the pipe. The slowing rate of the 
pressure pulse is shown by the increasing time duration of the pulse as 
the air percent increases. 2) Small amounts of air have a large effect in 
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reducing the amplitude of the pressure pulse and in decreasing the pressure 
pulse speed. 

The photographs taken at transducer location one show the same effects 
as the photographs taken at transducer location two. As a result, the 
photographs taken at transducer location one are not included in this re- 
port. Fig. 6 and 7 show the photographs taken at transducer location two. 
The results taken from the photographs are tabulated in Table 1. 

In Fig. 8 through 11, the theoretical curves include no frictional 
losses. Because there are no frictional losses included in the theoretical 
curves, it would be expected that the experimental results should fall 
slightly below the theoretical curve. Thi* is because the frictional 
losses for a two phase mixture are greater f\an those for a single phase. 
Martinelli experimentally showed this. / 1 0/ As discussed in the theory 

section, the compression of the air by the pressure pulse was probably best 
described as adiabatic because the method of injection allowed for large 
bubbles which gave a small heat transfer area per unit volume. The pres- 
sure pulse travels so fast that the air is compressed for a very short 
time, consequently this permits very little time for heat transfer. 

For the runs when no air was injected into the water, values of 81.7, 
71.7, 74.7 and 79.9 percent of the theoretical maximum pressure were obtain- 
ed at transducer location two for nominal heads of two, four, six, and 
eight feet respectively. Entrained air in the water of about 0.02 percent 
by volume would cause this effect. The entrained air was under pressure 
of about 100 psia in the United States Naval Postgraduate School system. 

When the water pressure was reduced, the air expanded and the volume per- 
cent increased. A glass of water will have air bubbles on the glass sur- 
face if allowed to sit overnight. 
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By allowing the water to stand about fifteen hours in the tank, the 
maximum pressure obtained when no air was injected increased from approxi- 
mately 50 percent to about 80 percent of the theoretical maximum. This 
indicates that the runs in which no air was injected had enough entrained 
air in the water that the exact air content was uncertain. A factor of 
0.02 percent by volume makes a large difference in the maximum pressure 
obtained in the region approaching zero percent of air. This same per- 
cent of uncertainty makes very little difference in the maximum pressure 
in the region greater than about 0.5 percent of air. Therefore, this 
entrained air did not effect the results appreciably for values greater 
than 0.5 percent of air. 

The curves of Fig. 8 through 12 were plotted on a scale in which the 

low ranges of pressure and volume percent were expanded. This was done to 

better show the region of greatest change. A log scale was used in which 

for the ordinate In ( P /P Nr + 0.1) was plotted and for the abscissa 

In (n + 1 .0) was plotted. This was done so that it would be possible to 

have zeros for P/P^, anc * ^ cn the scales. The results indicate that the 

points fall close to the theoretical curves in the regions greater than two 

percent of air. In the region of less than two percent of air, the results 

obtained are more erratic. The uncertainty in the region of 0.5 percent 

by volume air was plus or minus about ten percent in the volume percent 

term, and about plus or minus 11.0 percent in the P/P v _ term. Even consider 

NT 

ing these uncertainties, several of the points in the region less than two 
percent by volume air do not fall on the theoretical curve. In the region 
of greater than two percent by volume air the consistency of the data vas 
much better. At four percent by volume air, the uncertainty of the volume 
percent term (Y[ ) was about plus or minus 0.8 percent and that of the P/P^ 
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term was about plus or minus 17.0 percent. These uncertainties are 
based upon the uncertainty analysis of runs 3-9-2 and 3-9-10 given in 
Appendix H. 

Fig. 12 shows the effect of air in damping pressure pulses in systems 
at pressure up to 500 psia. This figure shows curves of 
for pressures of 500, 400, 300, 200, 100, and 15.5 psia. The volume per- 
cent of air at inlet conditions varies from zero to ten percent. These 
curves show that the lower the pressure in which the pressure pulse is 
being transmitted, the greater is the effect of air in reducing the pressure 
amplitude. As an example, at one percent by volume of air at 15.5 psia, 
the pressure would be reduced to 8.6 percent of the original pressure. 

At 100, 200, 300, 400, and 500 psi the pressure amplitude reduction would 
be to 20.8, 28.7, 34.4, 38.9, and 42.5 percent respectively of the original 
pressure amplitude. 

As an example of the usefulness of air in reducing pressure pulses 
and also to get an idea of the amount of air that would be required, con- 
sider the following example. 

Example 1. Pure water is flowing through a two inch extra strong 
pipe at the rate of 500 gpm. The system has an operating pressure 
of 100 psia and an operating temperature cf 70°F. Fluid borne noises 
are being transmitted in this system. It is desired to reduce the 
amplitude of these pressure pulses to 20 percent of the original 
amplitude . 

Fig. 12 shows that for the pressure amplitude to be 20 percent 
of the original amplitude, a 1.05 percent by volume mixture must be 
made by the injection of air. 
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.500,9 prfl- 
■7.48 gal/ffi 



— 46-9 4 ft /min water 



Let x equal the air required in ft /min at 100 psia and 70°F. 
Then 

0.0105 = 



100X 



X + 66.84 

x = 0.7093 ft 3 /min at 100 psia and 70°F 

3 

This equals 4.825 ft /min at 14.7 psia and 70°F. 

The weight of air required would be 0.3628 lbm/min. 
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7 . 



Conclusions , 



Based upon the experimental work and development of equations (6) 
and (6a) the following conclusions have been made. 

1. Small amounts of air in a water system greatly reduce the 
amplitude of a pressure pulse transmitted in the water. In a 15.5 psia 
water system with one percent by volume air, the pressure is reduced to 
8.6 percent of the theoretical no-air amplitude. 

2. The effect of air in damping pressure pulses decreases as 
the pressure of the system increases as demonstrated by Fig. 12. 

3. There was a reasonable correlation of experimental results 
with theory as shown in Fig. 8 through 11. 
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Fig. 1 System Diagram 



KEY TO SYSTEM DIAGRAM 



1. Water storage tank 

2. Two inch water fill line 

3. Storage tank overflow piping system 

4. Overflow drain tank 

5. Air injection system 

6. Two inch pipe coupling 

7. Fischer and Porter rotometer 

8. Potentiometer reference ice bath 

9. Air manifold 

10. Inlet air mercury manometer 

11. Thermocouple potentiometer 

12. Kistler Model 401 piezoelectric pressure transducer 

13. Kistler #470 shielded cable 

14. Kistler Model #651 piezo-calibrator 

15. Shielded calibrator-oscilloscope cable 

16. Hewlett Packard Model 130 oscilloscope 

17. Polaroid oscilloscope camera 

18. Falling weight guide tube 

19. Seven pound weight 

20. Modified Crane two inch quick closing valve 

21. Quick closing valve stop 

22. Water discharge mercury thermometer 

23. Dayton scales 

24. Weigh tank 
Air inlet themocouple 



25 . 
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26. Water inlet Weston thermometer 

27. Standard ASA two inch extra strong steel pipe 

28. One quarter inch vinyl tubing 

29. Two inch gate valves 

30. Air regulating valve 

31. Quick closing air stop valve 

32. Air flow regulating valve 

33. Water inlet regulating valve 

34. Weigh tank two inch discharge gate valves 

35. Oscilloscope sweep trigger 
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FIG. 2. Photograph of Air Injection System 





t?io. Fhotogranh of Valve and Closure Svstem 




KEY TO EQUIPMENT PHOTOGRAPHS 



Item 



1 . 



2 . 



3. 



A. 



5. 



6 . 



7. 



8 . 



9. 



10 . 



11 . 

12 . 



13. 



14. 



15. 



Description 
Thermocouple potentiometer 
Fisher and Porter rotometer 
Potentiometer reference ice bath 
Inlet air mercury manometer 
Air injection system 

Hewlett Packard model 130 oscilloscope 
Kistler model #6 5 1 piezo-calibrator 
Kistler #470 shielded cable 

Kistler model 401 piezo-electric pressure transducer 

Modified Crane quick closing valve 

Falling weight guide tube with weight in it 

Weigh tank 

Dayton scales 

Weigh tank discharge valve 
Electric timer 



24 




1 

,1 


I( 


I 1 

II 

1. 


II 

1 I 

L 




Fig. 5 Cross Sectional View of Air Injection Annular Ring 
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FIG. 6 PRESS'*!^ TIME 'HOTOGRAFHS PAPT 1 














6 

o 

c 

r: 



o 

r*? 

tj 

; < 

o 

tj 



I 

c 

o 



o 

-p 

o 



r. 

o 

C' 

•rt 

■rl 



O 



£ 



s 



CO 

E 



a 



*-< 

i 

fcU 

(h 



W 



co 

lu 



ffi 



t^- 

e> 

*— « 
t*- 
















28 



•J 

t; 




1.0 




d/d 





VO 




32 



Table 1. Tabulated Results 



Run No. 


Trans- 

ducer 

Location 


Head 

ft. 


Yl 

7. 


Y 

' WT 
7oxl0 


P 

psi 


P/P_ 

NT 


3-21-32 


1 


2 


1.09 


1.42 


10.6 


0.0498 


3-21-33 


1 


2 


1.09 


1.42 


10.6 


0.0498 


3-21-34 


1 


2 


3.05 


4.08 


10.6 


0.0498 


3-21-35 


1 


2 


3.10 


4.14 


11.2 


0.0526 


3-21-36 


1 


2 


5.14 


7.02 


8.8 


0.0413 


3-21-37 


1 


2 


5.17 


7.06 


8.2 


0.0385 


3-21-38 


1 


2 


7.19 


10.05 


7.3 


0.0343 


3-21-39 


1 


2 


7.19 


10.04 


7.3 


0.0343 


3-21-40 


1 


2 


9.23 


13.16 


6.4 


0.0300 


3-21-41 


1 


2 


9.28 


13.23 


6 . 4 


0.0300 


3-21-18 


1 


4 


0.664 


1.83 


22.4 


0.0704 


3-21-19 


1 


4 


0.664 


1.83 


22.4 


0.0704 


3-21-20 


1 


4 


2.02 


2.83 


14.6 


0.0462 


3-21-21 


1 


4 


2.01 


2.80 


15.8 


0.0497 


3-21-22 


1 


4 


3.37 


4.78 


11.1 


0.0349 


3-21-23 


1 


4 


3.37 


4.77 


11.9 


0.0374 


3-21-25 


1 


4 


4.77 


6.85 


10.4 


0.0327 


3-21-26 


1 


4 


6.38 


9.31 


10.4 


0.0327 


3-21-27 


1 


4 


6.38 


9.31 


10.4 


0.0327 


3-21-3 


1 


6 


0.538 


0.774 


30.6 


0.0765 


3-21-4 


1 


6 


0.536 


0.769 


28.8 


0.0720 


3-21-5 


1 


6 


1.53 


2.23 


26.4 


0.0660 


3-21-6 


1 


6 


1.53 


2.23 


24.6 


0.0615 
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Table 1. Tabulated Results (continued) 



Run No. 

: 


Trans- 

ducer 

Location 


Head 

ft. 


n 

% 


W WT 

%xl0" 3 


P 


P/Pxw. 

NT 


3-21-7 


1 


6 


2.60 


3.84 


18.0 


0.0450 


3-21-8 


1 


6 


2.60 


3.84 


18.0 


0.0450 


3-21-13 


1 


6 


3.68 


5.49 


15.6 


0.0390 


3-21-14 


1 


6 


3.68 


5.49 


15.6 


0.0390 


3-21-15 


1 


6 


4.76 


7.18 


12.0 


0.0300 


3-21-16 


1 


6 


4.80 


8.07 


14.4 


0.0360 


3-13-3 


1 


8 


0.475 


0.731 


35.8 


0.0773 


3-13-4 


1 


8 


0.475 


0.731 


41.7 


0.0901 


3-13-5 


1 


8 


1.28 


2.04 


33.5 


0.0724 


3-13-6 


1 


8 


1.28 


2.01 


34.6 


0.0747 


3-13-7 


1 


8 


2.14 


3.35 


25.7 


0.0555 


3-13-8 


1 


8 


2.14 


3.35 


29.3 


0.0633 


3-13-9 


1 


8 


3.01 


4.77 


19.1 


0.0413 


3-13-10 


1 


8 


3.01 


4.77 


20.9 


0.0451 


3-13-11 


1 


8 


3.91 


6.22 


17.9 


0.0387 


3-13-12 


1 


8 


3.91 


6.22 


19.1 


0.0413 


2-27-4 


2 


2 


0 


0 


174.0 


0.817 


4-1-4 


2 


2 


0.294 


0.382 


43.1 


0.202 


4-1-5 


2 


2 


0.294 


0.382 


45.4 


0.213 


2-28-4 


2 


2 


1.20 


1.65 


8.68 


0.0410 


2-28-5 


2 


2 


1.19 


1.58 


7.78 


0.0365 


2-28-13 


2 


2 


3.13 


4.25 


8.07 


0.0379 


2-28-14 


2 


2 


3.10 


4.22 


9.97 


0.0468 
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Table 1 Tabulated Results (continued) 



Run No. 


Trans- 

ducer 

Location 


Head 

ft. 


n 

% 


r/ 

' WT - 

7»xlU 


P 


P/P^ 

NT 


2-28-15 


2 


2 


5.20 


7.23 


5.67 


0.0266 


2-28-16 


2 


2 


5.18 


7.22 


6.27 


0.0294 


2-28-17 


2 


2 


7.20 


10.2 


6.03 


0.0283 


2-28-18 


2 


2 


7.28 


10.1 


5.37 


0.0252 


2-28-19 


2 


2 


9.54 


13.9 


5.97 


0.0280 


3-27-3 


2 


4 


0 


0 


228 


0.717 


4-1-6 


2 


4 


0.277 


0.257 


105.3 


0.331 


4-1-7 


2 


4 


0.277 


0.257 


84.1 


0.261 


2-28-26 


2 


4 


0.769 


1.11 


37.9 


0.118 


2-28-28 


2 


4 


2.087 


3.12 


19.9 


0.0626 


2-28-29 


2 


4 


2.087 


3.12 


15.7 


0.0494 


2-28-30 


2 


4 


3.49 


5.18 


12,1 


0.0380 


2-28-31 


2 


4 


3.51 


5.20 


12.7 


0.0399 


2-28-32 


2 


4 


5.02 


7.55 


9.10 


0.0286 


2-28-33 


2 


4 


4.99 


7.48 


12.1 


0.0380 


2-28-34 


2 


4 


6.24 


9.51 


9.1 


0.0286 


2-28-25 


2 


4 


6.20 


9.67 


10.9 


0.0343 


3-27-2 


2 


6 


0 


0 


228.1 


0.717 


4-1-8 


2 


6 


0.136 


0,207 


161.8 


0,404 


3-3-4 


2 


6 


0.596 


0.877 


28.6 


0.0715 


3-3-5 


2 


6 


0.602 


0.893 


38,8 


0,0956 


3-3-6 


2 


6 


1.57 


2.33 


16.9 


0.0422 


3-3-7 


2 


6 


1.59 


2.41 


26.2 


0.0655 


3-3-8 


2 


6 


2.59 


3.90 


22.0 


0.0550 
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Table 1 Tabulated Results (continued) 



Run No. 


Trans- 

ducer 

Location 


Head 

ft. 


n 

% 


3-3-9 


2 


6 


2. 


3-3-10 


2 


6 


3. 


3-3-11 


2 


6 


3. 


3-3-12 


2 


6 


4. 


3-3-13 


2 


6 


4. 


3-27-1 


2 


8 


0 


4-1-10 


2 


8 


0 . 


3-9-2 


2 


8 


0 . 


3-9-3 


2 


8 


0 . 


3-9-5 


2 


8 


1 . 


3-9-6 


2 


8 


2. 


3-9-7 


2 


8 


2. 


3-9-8 


2 


8 


3. 


3 9-9 


2 


8 


3. 


3-9-10 


2 


8 


3. 


3-9-11 


2 


8 


3. 



7x10 


P 


P/P*™ 

NT 


3.89 


16.0 


0.0400 


5.51 


17.8 


0.0445 


5.51 


17.8 


0.0445 


7.24 


15.4 


0.0385 


7.32 


10.6 


0.0265 


0 


370 


0.799 


0.179 


356.9 


0.771 


0.750 


51.9 


0.112 


0.747 


42.3 


0.0913 


2.12 


48.9 


0.106 


3.41 


26.7 


0.0576 


3.45 


27.9 


0.0602 


4.76 


19.5 


0.0421 


4.76 


21.3 


0.0460 


6.37 


17.1 


0.0462 


6.37 


15.5 


0.0462 



58 

61 

61 

70 

76 

117 

485 

484 

26 

13 

15 

01 

01 

89 

89 
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Fig. 13 Pipe and Tank Dimensions. 



APPENDIX B 



THEORY DERIVATIONS 

1, The following Is the derivation for equation (1) of the theory 
section as given by Stepanoff /13/. Newton's second law is as follows: 

tr — d(AAV) ( B - i 

dt; 

Where: F = Force, lbf. 

V = Velocity, ft/sec. 
t = Time, sec. 

M = Mass, slugs. 



* — X- <3 t • 




The velocity through the pipe (V ) before valve closure is constant. The 
force (increase in pressure at the valve) is equal to the time rate change 
of momentum of the moving mass. The pressure- force relationship is as 
follows : 

F = A v H 3 (b-2) 

At time t. after closure the portion of the fluid in the length of pipe, X, 
will have its velocity reduced to zero, then: 

F=A<TH,= V. <»- 3 > 

2 

Where: A = Pipe internal cross sectional area, ft 

Y = Specific weight of fluid flowing, lbm/ft^ t 
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H = Head rise above initial no flow static head, H , ft of 
y o 

J water 

g = Proportionality constant in Newton's second law of motion, 
(lbm-ft)/ (lbf-sec2) 

X = Distance to any position on the pipe measured from the 
valve, ft. 

t^ = Time after closure, sec. 

V q = Velocity of fluid in the pipe before valve closure, ft/sec. 

The term A — in equation (B-3) is the mass flow rate. In equation 
“o 

(B-3), X/ 1 ^ is the velocity of the pressure pulse, a. Substituting a and 
cancelling the Y A term on both sides of the equation (B-3) gives: 



Li ^ Vo 

= a -57 (1) 

2. The derivation of equation (2) in the theory section as given by 
Stepanoff / 1 3/ , is as follows. After instantaneous valve closure a region 
of compressed water is created adjacent to the valve and extends a distance 
X at time t^ after closure. The pipe in this region expands under pres- 
sure from D to D- . From the instant of closure until t. a volume of wa- 
o 1 i 

ter equal to AV^t^ passes through the original cross section at X. There 
is no outlet for the water; therefore, the volume must be accommodated by 
expanding the pipe radially and by compressing the water. 




hocus 
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The increase in volume caused by the compression of the water in the pipe 
is obtained from the definition of bulk modulus as follows: 



Kw cHG 



dP _ r 

i T7T O w jl ^ 



W 



AG 



(B-4) 



w 



Where: K = Water bulk modulus, lbf/ft . 

w 

3 

= Volume, ft 

AGy = Incremental change in volume, ft' 



P = Incremental change in pressure caused by the velocity head 
changing to pressure head, lbf/ft • 



Solving equation (B-4) for AG W gives: 

= ~ G * £- 
yv 



(B-5) 



but the original volume G„ is equal to AX. Substituting AX for G„ gives: 

w W 



AGw = -AXJL 

* w 



(B-6) 



The negative sign indicates that AG decreases with increasing pressure. 
The decreased volume is used to accommodate the additional mass that 
enters the section of Pipe X feet long. The volume made available by the 
compression of the original water is: 

_ A X P 



AG W - 



K 



(B-7) 



w 



The increased volume caused by the expanding of the pipe is: 



AGp = f(D. i -D/JX 



(B-8) 



By equating the additional volume of water that must be accommodated 
(AV^t^) to the volume that the pipe expands plus the volume that the 
water compresses the following is obtained: 
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Av.ti - -45 + ^(o;--o/jx 



For the pipe 




PO e 



P D 0 = 2Se 



Where: 

e = Pipe wall thickness, ft, 
Hooke's law for the pipe expansion is: 

D.-O, 



-E( 



Dc 



Where: E = Modulus of elasticity, lbf/ft . 

Combining equations (B-10) and (B-ll) gives: 



(B-9) 



(B-10) 



(B-ll) 



P Do 
2 aE 



P,-Oo 

Do 



( B - 1 2 ) 



2 2 

In equation (B-9) -D may be written (D^+D q ) (D^-D q ) . 

Since Dj' D 0 is very small compared to or D q , the term (D^+D^) may be 

written 2D q . Then the last term of equation (B-9) becomes: 

_ ~x.fr P XP.-0-) (B . 13> 

4 ~ z 

Multiplying both sides of equation (B-12) by gives: 



d ( -o. = JLSfL 

2eE 



(B - 1 4 ) 



Substituting equation (B-14) into equation (B-13) and substituting A for 

gives 

x it ( p Dc 1 ) _ 2L a PqP 
4- eE 



the cross sectional area term, 



(B- 15) 



Substituting equation (B-15) for the last term in equation (B-9) gives: 
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V. A t; 
V.t i 



X PA 

XP(± 



H- 



x Oo pa 

e. E 

Qo \ 

eEJ 



(B-16) 



The term X/t^ is the wave velocity, a. Substituting and solving for a 



gives : 



Q = 



Vo 



( -L 4- JBi) 
V K w eE/ 



(B- 17) 



If equation (1) is multiplied by the specific weight of water ( "$" ) , the 



following is obtained: 



1 1 v o ^ v» 

2T W — r — cj 0 



(B-18) 



Solving (B-18) for V q gives: 



V. = 



Pc 



Q r 



w 



Substituting (B- 19) into (B -17) and letting D equal and J ^ equal 
gives : 



(B-19) 

X, 



a = 



+• ek) 



( 2 ) 



The maximum theoretical pressure rise, (P^p) * s given by substituting equa- 



tion (2) into equation (B-18). 

P.X = 



V. V, 






tTW) 



(4) 



The experimental results were expressed as P/ P^, 



P _ P '/p-C^L * ~fe) 






Vo 



(B-20) 



The term P is the experimental pressure rise above the no flow static pres- 
sure * 
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3. The following is the derivation of equations (3), (3a), (5), (5a), 

(6), and (6a) of the theory section. These derivations are based upon 
work done by Karplus /4/. His equations were modified to make them apply 
to the water hammer model. The following is a form of equations (2) for 
an air-water mixture. 




(B-21) 



Where: Subscript c represents the air-water mixture. 




c 



2 

mixture compressibility, ft /lbf. 



The density of the mixture ( 




c 



) is given by: 




(B-22) 



3 

Where: G = Volume^ ft . 



M = Mass, slugs. 



subscripts 



a = Air 



W = Water 



Because mass is the product of density and volume ,0 becomes : 




(B-23) 



The volume fraction of air is defined as: 




(B-24) 



The volume fraction water is defined as: 



^ = -§? 



(B-25) 



Substituting (B-24) and (B-25) into (B-23) gives; 




(B-26) 
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Compressibility is defined as follows: 



c\j _ d£c. —(olCg 0I6J 

Ac - G c dP ~ G c dP 

d G vy /Gw/ \ dG a* ( Ga \ 

~ S w dP UJ \gJ 

Substituting equations (B-24) and (B-25) into (B-27) gives: 

X = 0-*) X, + x ^«w 

Substituting equations (B-26) and (B-28) into equation (B-21) gives: 



Rearranging equation (B-29) gives: 



Q' l = +( 1 -> 0 V’wX+ X(l-x)f/.X+/^^n]+£ E §-«)/rt ♦*/>,] 



For isothermal conditions and assuming the ideal gas law applies: 



i_ 

Pass 

Where: ^ABS = Absolute Pressure of system, lbf/ft 

^ = Isothermal compressibility of air. 

The velocity of sound in air is given as: 




2 

ft 2 /lbf 



Where %a= Adiabatic compressibility of aii^ ft /lbf. 

For adiabatic conditions and assuming the idea'l gas law applies: 

__ i 

A Aft ft, p 

’A 8S 



(B-27) 



(B-28) 



(B-29) 



(B-30) 



(B-31) 



(B-32) 



(B-33) 
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Where : 



k = ratio of specific heat of air at constant pressure to that 



at constant volume, 1.403 cor air 

By combining equations (B-31) and (B-33), the following is obtained: 

<»- 34 > 

The velocity of sound in water is given by: 

= /w “K*, 

For isothermal conditions substitute (B-32), (B-33), (B-34), and (B-35) into 
(B-30). This gives: 



Q 



-2 
X ~ 





+■ 






(3) 



Using the adiabatic conditions for equations (B-30) gives the following 
results : 




(3a) 



For isothermal conditions for an air-water mixture equation (B-18) becomes: 




Vo Qr 



(B-36) 



Where: a^ is given by equation (3). 

Substituting equation (B-36) into (B-20) gives: 




( 6 ) 



Where : a^. is given by equation (3). 
For adiabatic conditions: 
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(6a) 




Where: a is given by equation (3a). 

A 
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APPENDIX C 



EQUIPMENT 

1. The head tank consisted of a tank which was two feet in diameter at the 
lower end and one foot ten and one-half inches at the upper end. The upper 
part of the tank was light sheet metal and the lower part of the tank was 
plate steel. Two inch overflow lines were provided at approximately two 
foot intervals as is shown in Appendix A. These overflow lines were stan- 
dard two inch galvanized steel pipe. Each overboard line had a two inch 
gate valve in the line. All the overboard lines ran into a common line which 
discharged into a collecting tank open to the atmosphere. This collecting 
tank had a drain which discharged into the main sewage system, A Weston 
thermometer located near the bottom of the tank, was used to measure the 
inlet water temperature. 

The main run of pipe consisted of two inch extra strong black steel 

pipe. The out side diameter was 2,375 inches, the inside diameter was 1.939 

inches. The wall thickness was 0.218 inches. Three sections of pipe were 

welded together. Care was taken to insure good butting of the individual 

♦ * 

sections. At the transducer locations, transducer fittings one inch in 
diameter and one inch high were welded into the pipe. These fittings had 
standard 14 mm by 1,25 automotive spark plug threads. 

The water we*gh system consisted of a Toledo Scales Company platform 
scales with a capacity of 1600 pounds and a steel plate rectangular tank, 

20 inches by 28 inches by 20 inches. A 3 inch by 3 inch section was cut 
out of the top to permit entrance of the discharge pipe. This tank had 
two, 2 inch gate valves for discharge. They discharged into drain lines to 
the sewage system. 
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The following scale information is provided. 



Manufacturer 



Toledo Scales Company 
Toledo, Ohio 



Model 



31-1811FC 



Serial Number 



2180-0-225 



Capacity 



0-1600 pounds 



An electric timer, which measured to the nearest one tenth cf a 
second, was used to measure the time to fill the weigh tank with water. 

2. The air injection system is shown schematically as part of Fig. 1. 

A photograph of the air injection system is shown in Fig. 2. The source of 
air was from the 125 psi air system permanently installed in Building 500, 
at the United States Naval Postgraduate School. This system was reduced 
to approximately fifteen psi by a regulating valve. One-quarter inch vinyl 
tubing was used from the regulating valve to a pipe manifold. A copper- 
constantan thermocouple was installed in the pipe manifold to measure the 
temperature of the air. A Rubicon Company catalog number 2732 potentiometer 
was used to measure the thermocouple output. This is item 1 in Fig. 2. 

To measure the pressure of the air a 50 inch mercury manometer was used. 

A one-quarter inch vinyl tubing line was tapped into the air manifold. This 
line was connected to the mercury manometer. The other side of the manometer 
was open to the atmosphere. Item 4 of Fig. 2 shows this manometer. From 
the manifold a one-quarter inch vinyl line led to the rotometer. The follow- 
ing data is provided. 

Manufacturer Fisher-Porter Company 

Flowrotor tube number 2-F-1/4-20-5/70 

Float glass 



Range (14.7 psi, 70°F) 



0-0.325 ft 3 /min. 
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The manufacturer's calibration curve, converted from cubic centimeters per 
minute to cubic feet per minute, was used. Fig. 14 is a copy of this 
curve. Four points on this curve were checked using American Meter Com- 
pany wet gas meter Al-20-11472. 

The flow rate was corrected to 14.7 psia and 70°R for each run using 
the following relationship suggested by the manufacturer / 2/ . 

Q = ft'Vmin at 70°F, 14.7 psia 

3 

= pressure at which measurement was made, psia 

= temperature at which measurement was made, °R. 

3 

Q = flow rate in ft /min taken from the curve of Fig. 14. 

This curve was entered with the rotometer reading. 

Item number 2 in Fig. 2 shows the rotometer. A one-quarter inch vinyl 
tubing connected the rotometer with a manual control valve. The control 
valve was used to control the air flow rate. From the control valve a 
one-quarter inch vinyl tube led to ; manifold section in the pipe. A 
stop valve was provided for each of these separate lines. Item 5 of Fig. 

j 

2 shows the equipment used for the injection of the air into the pipe. 

Three separate rings of holes were used to inject the air. One ring con- 
sisted of two holes, each with a diameter of 0.04 inches. The second ring 
consisted of four holes, each with a diameter of 0.07 inches. The third 
ring consisted of eight holes with diameters of 0.07 inches. These holes 
were equally spaced around the circumference of the pipe. A plexiglass 
ring, four inches in diameter, with a circumferential groove for air injec- 
tion was provided for each ring of holes. These rings were sealed with "0" 
rings. Fig. 5 shows a cross sectional view of this system. The system was 
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designed so that any one of the rings could be used. However, it was 
not necessary to do this. All experimental work was done with air enter- 
ing all three rings of holes. 

3. The valve and closure system is shown in Fig. 3. The quick closing 
valve was a modified standard Crane two inch quick closing gate valve. 

The valve was modified by screwing on a steel disc two inches in diameter 
and one inch high on the top of the stem. This disc was the point of im- 
pact for the falling weight which closed the valve. A mechanism was de- 
signed to hold the valve closed when it was in the shut position. A stand 
that was completely isolated from the piping system was used to stop the 
valve. By using this method, the full impact of the stopping valve was 
not transmitted to the piping system. A seven pound cylindrical steel 
bar one inch in diameter by 24 inchs long was used as the weight for clos- 
ing the valve. This weight was dropped through an aluminum pipe, one 

and one-quarter inches in diameter. This aluminum pipe contained numer- 
ous air holes to allow the weight to fall freely. Appendix G contains pre- 

liminary calculations used to estimate the closure time. The weight was 
dropped from nine feet for all runs. The closure time of the valve was 
measured with an oscilloscope. The time sweep of the oscilloscope was 
triggered by a photoelectric cell. A contact completed a 45 volt direct 
current circuit when the valve was closed. This voltage was put across 
the vertical plates of the oscilloscope. A Polaroid time exposure was 
made of the oscilloscope, the closure time was taken from the photograph. 
Polaroid 3000 speed film was used. 

The quick release mechanism consisted of a steel pin which held the 
weight nine feet above the valve. A line was attached to this pin. When 
it was desired to drop the weight, the pin was pulled and the weight fell. 
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A nylon string was attached to the weight. This was used with a 



pulley to pull the weight to the desired height above the valve. 

4. Selected characteristics of the pressure sensing equipment are given 
below. For complete details refer to manufacturer’s instruction book / 5/ . 

Kistler 401 Pressure Transducer 



pressure range 


0.0-3000 psi 


overload capacity 


507. 


linearity 


17. 


natural frequency 


48 s 000 cps 


maximum continuous 
temperature 


500°F 


maximum intermittent 
temperature 


3000 °F 


length 


2.04 in. 


diameter 


0.625 in. 


weight 


1.75 oz. 


serial number 


#1044 



Kistler Model 651 Piezo Calibrator 



input voltage , maximum 


0.6 volts 


output voltage, maximum 


0. 3 volts 


voltage gain 


0.5 


input impedance 


,„14 

10 ohms 


output impedance 


lO^ohms 


frequency range 


0-10 kc 


linearity 


17. 


dimensions 


7 x 10 x 6 in. 


weight 


4.7 lbm 
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Hewlett Packard Model 130 Oscilloscope 



internal sweep accuracy 5% 

Vertical and Horizontal Amplifiers 
frequency response 0-300 kc 

sensitivity lmv/cm-20 volts/cm 

accuracy 57 0 

input impedance 1 megohm 

Cathode Ray Tube 

power requirements 115/230 volts, AC 

50/400 cps, approximately 
175 watts 

dimensions 9 3/4 in. by 15 in. by 21-1/4 in. 

weight 39 lbm 

Kistler transducer model 401 serial number 1044 was calibrated using 
an Ashcroft dead-weight tester* The system used for calibration was the 
same as shown in Fig. 4 except the transducer was inserted into the 14 mm 
dead-weight adapter fitting instead of the pipe. Sweeney / 1 4/ describes 
the procedure for using the dead-weight tester. The important thing in 
calibrating was to rotate the weights- when loading the transducer. The 
transducer was rapidly loaded so that there would be no error due to 
charge leakage from the crystal. The transducer was calibrated for the 
calibrator range scales of 10 and lOf.The 10 range covers 0-125 psig and 
100 range covers 0-1250 psig. Fig. 15 and 16 show these calibration 
curves. These curves show the readings on the calibrator for a given 
pressure. It was necessary to determine the number of calitration divisions 
that correspond to a centimeter deflection on the oscilloscope. Knowing 
this and using the calibration curves* the psig/cent imeter on the oscillo- 
scope was calculated. For the oscilloscope used for this experimental 
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PIG. 15* Transducer Calibration Curve, 0 - 125 ps’ig 
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PIG# 16# Transducer Calibration Curve, Q-1250 psig 




work, the following values were obtained and used. These values are for 



Kistler model 401 serial number 1044 pressure transducer,, 



Calibrated 

range 


Osci llcseopa 
scale 


Calibration 

divisions/ 

10 cm oscillo- 
scope deflection 


psig/ cm 


10 


5 mv/ cm 


19.7 


3.0 


10 


10 mv/cm 


39.1 


6.0 


10 


20 mv/cm 


77.2 


11.8 


100 


5 mv/cm 


17.2 


26.6 


100 


10 mv/cm 


34.4 


54.4 


100 


20 mv/cm 


67.8 


106.9 
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APPENDIX D 



DETAILED PROCEDURE FOR A RUN 

The detailed step by step procedure for each run was as follows: 

1. Record atmospheric temperature and pressure. 

2. Fill the tank to the level desired. 

3. Open quick closing stop valve. 

4. Adjust water inlet to the tank until the overflow at the 
desired level is just a trickle. 

5. Adjust the air rotometer to desired flow rate with the 
^control valve immediately following the rotometer in the air 

inlet system. 

6. Insure weigh tank discharge valves are completely open and 
that the tank is as low as possible. 

7. Set electric timer to zero. 

8. Close weigh tank discharge valves. 

9. Start timer and record initial weight on scales when the 
timer is started. 

10. Fill weigh tank until it is nearly full. 

11. Record weight and time at nearly full position. 

12. Open weigh tank discharge valves. 

13. Lift weight to nine foot position and insert quick release 
pin. 

14. Record inlet and outlet water temperatures, 

15. Record rotometer setting (read ball top). 

16. Record air mercury manometer readings. 

17. Record air temperature. 
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18. Set calibrator range to desired position. Insure 

the meter is adjusted to 1,0. Allow a ten minute warm up 
period for the calibrator and the oscilloscope, 

19. Set oscilloscope on desired vertical scale. 

20. Adjust oscilloscope sweep rate to the desired time scale. 

21. Insure the SYNC switch is set to external trigger (EXT) on 
the oscilloscope. 

22. Insure that battery is properly connected to trigger 
circuit of the oscilloscope. 

23. Set the horizontal sensitivity of the oscilloscope to 
internal sweep. 

24. Check system by checking the calibration reading for ten 
centimeters deflection on the vertical scale of the oscillo- 
scope. 

25. Adjust vertical and horizontal position until the oscillo- 
scope trace is in desired position. 

26. Set camera on time exposure. 

t ¥ 

27. Open camera shutter. 

28. Close valve by dropping the weight with the quick release 
mechanism. 

29. Close camera shutter. 

30. Open quick closing valve. 

31. Develop Polaroid photograph. 
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HORIZ NTAl. EXPERIMENTAL OATA SHEET 




APPENDIX F 



SAMPLE CALCULATIONS 

Sample calculations have been included for run number 3-9-10. The 
data for this run has been included on the sample data sheet in Appendix 
E. 

P ATM = 30 • 264 in * H S 

t«™ = 76. 0°F 
ATM 

c = -0.129 Barometer temperature correction 
Calculate atmospheric pressure in psia 

( p ATM ( in «8) + c) (0.4913) = p ATM (P sia ) 

(30. 264-0. 129)(0. 4913) = 14.81 psia 

where 0.4913 is the conversion factor for converting pressure in inches 
of Hg to psia. 

Calculate absolute pressure of air. 



(P l 3 + ? 2a ) (0 ' 4913) + P ATM = p a (P sia > 

P^ a and are the inches of mercury measured on each leg of the air 

manometer. 

P la = 12 ’ 62 in * Hg 
P 2a = 12.55 in. Hg 

P ATM * U ' 81 pSi “ 

= pressure at which air was measured, psia 
(12.62 + 12.55)(0.4913) 4- 14.81 = 27.18 psia 
Calculate air flow rate at 14.7 psia and 70°F. 
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where: 



Q = air flow rate at 14.7 psia and 70°F, ft /min. 

d 

Q = uncorrected air flow rate taken from Fig, 14, 
u 

Enter Fig, 14 with rotometer setting. 

= Absolute temperature at which air flow rate was measuredj 



° R / 

/ Sio \ 

Q a -- °- 325 J LtjzJ 



= 0.442 ft^/min 



Calculate mass flow rate of air 



™a = Q a X* 



where : 

= Mass flow rate of air, lbm/min 

y = Specific gravity of air at 14.7 psia and 70°F. 
a 

ma = (0.442) (0.0752) = 0.03403 lbm/min 
Calculate volume flow rate at inlet water conditions 



= Q a( J r£) (yk) 



where 



Q = Air volume flow rate at conditions of water inlet 
aw 

= Absolute pressure of water at air inlet, psia 
= Absolute temperature of air at water inlet, °R 



Q = 0.442 
aw 



(J±l )(. 

V 18.245/ L 



£i± ) 

S3oJ 



= 0.3463 ft /min 



Calculate the water mass flow rate 

Yrt .= foo) 

v-J 4. • 

• 

Yfl = Water mass flow rate, lbm/min 

• ' W 

VY1 W = Weight of water weighed in time t^, lbm 
t^ = time to measure pounds of water, sec 
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m tT = = 534 lbm/min 

w 21.4 

Calculate water volume flow rate 



where : 



*W 



ra w 



3 

Q._ = Water volume flow rate, ft /min 
W 

Y w = Specific weight of water, lbm/ft 

Q W = ^624“ = 8 - 56 ft3 / min 

Calculate the volume percent air in mixture 




Q 

Q 



aH < 100 > 

aU + % 



where : 



= Volume percent air at water inlet conditions % 

" ’ (0.3463 + 5.56) ~ 

Calculate the weight percent air in the mixture 

vm) 

Wl,^ +■ 



where : 



n 



WT 



= Weight percent air in air water mixture, % 



tf (0.03403) (100) 

U W7 ~ (0703403 + 534) = 0 * 006377 - 



Calculate the maximum pressure 
P = d#C-J 

P = Pressure above the static no flow pressure, psi 
C = Psi/cm of oscilloscope deflection, psi/cm 
d = Deflection of oscilloscope, cm 

J = Correction for velocity head at the location when the 
oscilloscope scale was zeroed, psi. 
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To determine the J correction, the difference in static pressure with flow 
and with no flow was experimentally measured for each head at each loca- 
tion. It was necessary to make this correction because the oscilloscope 
was zeroed for each run when there was velocity head in the pipe. 

P = (3. 1) (6.0)-2, 1 = 16.5 psi 
Calculate pressure ratio P/P^ 



P/P 



NT 



P/P 



.,/£7I7K) ’• 



= P ao . v ^ y( ' L .i&) 

NT w v 



where : 



Vo* w 



V q = Velocity of the water through the pipe with no air, 

determined experimentally, ft/sec 

2 

= Density of the water, slugs/ft 

2 

Bulk modulus of water, lb/ft 
D = Inside diameter of pipe, in. 
e = Pipe wall thickness, in. 

2 

E = Modulus of elasticity of pipe lbf/ft 

= Theoretical maximum head rise above the initial static 
head, equation l,ft of water 



) 



i.sis 



p/p = ( 3l0 /° oO Xl44) 4 (o.2i*)(E3.4» 109044) 

m (?.SZ)(02.4) 



P/P„^ = 0.0445 
NT 
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APPENDIX G 



VALVE CLOSURE 

The water hammer model of Fig. 1 was used to generate the pressure 
pulses. For the maximum pressure to be theoretically possible, the model 
required a valve which would close in less than — seconds, 

where: L = The length of the pipe from the tank to the valve, ft. 

X = The distance from a given location to the valve, ft. 

a = The velocity of the pressure pulse in the pipe, ft/sec. 

The speed of the pressure pulse in pure water in a two inch extra strong 

steel pipe is 4580 ft/sec. For the maximum pressure as given by equation 
(1) to be possible at transducer location two, the valve had to close to 
less than 0.0128 seconds. If there was 0.05 percent by volume air at the 
water inlet conditions tie closure time required at transducer location 
two for the maximum pressure to be possible was 0.0351 seconds. With ten 
percent by volume air, the closure time required at transducer location 
two was 0.438 seconds or less. Since the pressure wave travels fastest 
in water with no air, this was the limiting condition. Therefore, to use 
transducer location two, a closure time of less than 0,0128 seconds was 
required. 

A free falling seven pound weight starting with zero initial velocity 
and falling nine feet would have a velocity calculated as follows: 

v = y"2 g s 

V = v/2 (32 , 2) (2) = 24.0 ft/sec 

2 

where g = acceleration of gravity, ft/sec 
s = distance the weight falls, ft 
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V = velocity of the falling weight, ft/sec 
With a 24.0 ft/sec. velocity, the time of closure would be: 




where t^ = time for valve gate to travel S v distance, sec 
S v = distance valve gate travels for closure, ft 
V = velocity of valve gate 
t t = (1/ 6) (1/24) = 0.0069 sec. 

Assuming an plastic collision, the following calculations are made. 



M 1 V 1 = (M 2 + M 1 )V 2 
M 1 

V = — V 

2 M 1 +M 2 1 

Assuming M 2 = 1.75 lbm = 1.75/32.2 slugs 

Mj = 7 lbm = 7/32.2 Slugs 

/_7 1 (24) = 19.2 ft/sec. 

V 2 = \(7+l. 75)^ 

The closure time for the valve would be: 

fc i = S i /V 2 = (6) (19.2) = °* 0087 sec * 



The packing on the valve was loose enough that the valve closed freely 
with very little resistance. The weight of the valve stem and gate was 
1.75 lbm. 

Experimentally the closure time was measured as 0.007 seconds with an 

uncertainty of 0.003 seconds. The method of measuring the closure time 

♦ 

lias been discussed in Appendix C. 
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APPENDIX H 



UNCERTAINTY ANALYSIS 

The uncertainty for runs 3-9-2 and 3-9-10 are included here as re- 
presentative examples of the uncertainty of the experimental work. Since 
the magnitude of the variables differs between runs, these values are to 
be considered only as representative. The method of analysis is that pro- 
posed by Kline and McClintock / 6/ for single sample data. 



Run Number 


3-9-2 


% 


3-9-10 


7« 


Item, Units 


Value + Uncertainty 


Value + Uncertainty 


m , lb/min 

V 


563.6 + 13.6 


2.4 


534 + 12 


2.2 


3 

Q , ft /min 
w 


9.03 + 0.22 


2.4 


8.56 + 0.19 


2.2 


m , lb 
w* 


170+ 4 


2.4 


190 + 4 


2.1 • . 


water 
t., sec 


18.1 + 0.1 


0.6 


21.4 + 0.1 


0.5 


T , °R 
w 


516 + 1 


0.2 


518 + 1 


0.2 


P , psia 
w 


18.26 + 0.07 


0.4 


18.26 + 0.07 


0.4 


3 

Q , ft /min 
u 


0.040 + 0.004 


1.0 


0.325 + 0.004 


0.1 


3 

Q , ft /min 

cl 


0.0562 + 0.006 


10.0 


0.442 + 0.0057 


1.3 


3 

Q , ft /min 
aw 


0.04403 + 0.00439 


9.9 


0.3463 + 0.0047 


1.4 


P, psia 


51.9 ± 6 


10.6 


17.1 + 3 


17.5 


P/P NT' 


0.1121 + 0.0129 


10.6 


0.0368 + 0.0064 


17.5 


P atm> "S 


30.264 + 0.002 


0.007 


30.264 + 0.002 


0.007 


T atm’ ° F 


76.0 + 0.5 


0.65 


76.0 + 0.5 


0.65 


ft, * 


0.485 + 0.050 


10.2 


3.89 + 0.77 


0.8 
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